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Abstract 

The maximum absorption of solar radiation over the broadest range of frequencies and inci- 



, ^ dent angles using the thinnest material possible has important applications for renewable-energy 

^^ length of scattered light rays, first proposed nearly thirty years ago, has thus far remained the 

c/3 most effective approach for photon absorption over the widest set of conditions. Recent thin- film 
O 



nanostructured designs involving resonant wave effects of photons have explored the possibility of 

superior performance though as of yet no proposal satisfying the dual requirements of enhanced 

and robust absorption over a large fraction of the solar spectrum has been made. Here using recent 

^* advances in computational electrodynamics we describe a general strategy for the design of a sili- 



con thin film applicable to photovoltaic cells based on a quasi-resonant approach to light trapping 



psj where two partially-disordered photonic-crystal slabs, stacked vertically on top of each other, have 

> 

On large absorption that surpasses the Lambertian limit over a broad bandwidth and angular range. 

(N 

^n This tandem structure, the photonic analogue of the multi-junction cell, has almost double the 

Tjj- improvement in the light trapping relative to a single-lattice design with equivalent thickness and 

o 

en is five times more effective than an unpatterned slab. 
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Photovoltaic (PV) cells that can efficiently and cost-effectively convert incident solar 
radiation into electricity have the potential to become an important component of appli- 
cations for renewable-energy generation pp. One of the fundamental issues underlying the 
design of silicon cells for use in realistic settings is the maximum absorption of incident 
sunlight over the widest possible range of wavelengths, polarizations and angles using the 
thinnest material possible. While complete random or so-called Lambertian texturing of 
the surface to isotropically scatter incident light rays into a weakly-absorbing thick film so 
as to increase the optical path length, as shown in Fig. flla), has thus far remained the 
most effective approach for light trapping over a wideband spectrum [21 E], recent thin-ffim 
nanostructured designs including photonic crystals [1] shown in Fig. U[b) exploiting the 
more complicated wave effects of photons have explored the possibility of superior perfor- 
mance [514T7] but have been mainly limited to narrow bandwidths, a single polarization or a 
restricted angular cone typical of delicate resonance-based phenomenon. The introduction 
of strong disorder, while improving robustness, nevertheless comes at the expense of light 
trapping relative to the unperturbed case [ISl US]. As a result, no proposal derived from a 
conceptual framework of photon capture for a nanostructured silicon thin film capable of 
robust, super-Lambertian absorption over a large fraction of the solar spectrum has yet been 
made. In this work, we present a new approach to light trapping made possible by recent 
advances in computational electrodynamics |20] based on the quasi-resonant absorption of 
photons that combines the large absorption of optical resonances with the broadband and 
robust characteristics of disordered systems via a stacked arrangement of ordered PC slabs 
augmented with partial disorder that surpasses by a wide margin the performance of an 
idealized Lambertian scatterer over a broad spectrum and angular cone. 

Our tandem design, consisting of the same silicon film structured in two different ways and 
stacked vertically as shown in Fig. fltc), is the photonic analogue of the multi-junction cell 
that employs three or more different semiconductor ffims where the electronic bandgaps add 
complimentarily to obtain wideband absorption. Here we demonstrate the utility of a pho- 
tonic approach, employing geometric structure alone, to enhance light trapping which also 
offers improved performance but does not incur the constraints and limitations of optimally 
combining multiple material-specific non-tunable electronic bandgaps as well as the signif- 
icant fabrication challenges and costs of epitaxially growing films with mismatched atomic 
lattices on top of one another. We outline a two-part design strategy based on first maxi- 



mizing, with a few-parameter gradient-free topology optimization the number of resonant- 
absorption modes by using two crystaUine-sihcon PC slabs with a fixed total thickness of 
1/im stacked on top of each other such that their individual resonances add compliment arily 
over the wideband spectrum; and then in the final step introducing a partial amount of 
disorder to both lattices to maximize the overall light trapping and boost robustness to go 
well beyond the Lambertian limit. 

In our earlier work, we showed how individual resonant-absorption peaks of a thin-film PC 
slab can be broadened using partial disorder leading to an overall enhancement of the wide- 
band absorption spectra and based on this proposed the following heuristic design strategy: 
find the maximum amount of disorder that still preserves some fundamental Bloch character 
of the modes and retains features of the large peak amplitudes to improve light trapping and 
robustness without broadening too much so as to be completely incoherent and leaky [TS] . 
To understand quantitatively why disorder increases broadband light trapping in a PC, we 
use coupled-mode theory to derive an analytical expression for a single absorption resonance 
(at a frequency of uq) in terms of the decay lifetimes (proportional to the quality factor) for 
radiation (Trad) and absorption (Tabs) by the material [1]: 



A{uj) = ^^'i'^^^ 2- (1) 
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When the decay lifetimes are matched, the absorptivity is maximized to one at the 
resonant frequency which is useful for thermal devices where narrowband, large emissivity 
is desirable [21j. However, broadband absorption for a thin-film PC design, consisting of 
a collection of such individual Lorentzian peaks [221 ES] , necessitates that we consider the 
total area spanned by eq. ([T|: 

A{u)du = — . (2) 
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The effect of disorder, as we have already demonstrated |il8j , is to reduce the peak height 
but more importantly to broaden the peak width (proportional to l/Trad+^/Tabs) via primar- 
ily a decrease in Trad which therefore leads to an overall increase in broadband absorption 
from eq. ^ (though Tabs also changes with disorder due to variations in the nature of 
the guided mode, the change is much less pronounced than Trad mainly because the mate- 
rial absorption coefficient is fixed). Note that this analysis is only valid when coupling to 



a resonant Bloch mode which is why introducing too much disorder and ehminating the 
peaks altogether, thus transitioning to non-resonant Anderson-locahzed modes, resuhs in 
sub-optimal light trapping [IB]. 

We consider the absorption of solar radiation in the wavelength regime spanning 600nm 
to llOOnm in which silicon is poorly absorbing and thus weak coupling to resonant Bloch 
modes of the PC most apparent. The absorptivity of PC thin films below wavelengths of 
600nm is usually quite large and little affected by the choice of lattice, our principal design 
objective, and thus omitted from the analysis while beyond llOOnm silicon's absorption 
coefficient is insignificant. We quantify the overall light-trapping efficiency of each design 
relative to an ideal perfect absorber which has unity absorptivity over the wavelength interval 
of interest by assuming that each absorbed photon with energy greater than the bandgap of 
silicon generates an exciton which contributes directly to the short-circuit current (properly 
accounting for the effects of charge collection in the device is unnecessary here given that 
light trapping as it relates to the Lambertian limit is purely an optical phenomena). This 
corresponds to the following definition of photon-absorption efficiency: 

Cr«:(A)-4(A)dA 
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where X(A) is the terrestrial power density per unit wavelength from the sun at AMI. 5 [21] 
and ^(A) is the absorptivity of the film. 

The design strategy of maximizing the light-trapping efficiency by controllably introduc- 
ing a partial amount of disorder to obtain just the right dose of peak broadening ultimately 
results in a more-uniform absorption profile where the absorptivity in the inter-peak re- 
gions is increased at the expense of the height of all peaks. This therefore suggests that 
in order to most effectively make use of partial disorder for light-trapping enhancement in 
a thin ffim the number of resonant modes must ffist be made as large as possible so that 
there is little bandwidth separation between peaks: by extending our previous slab design 
to include not one but two PC lattices stacked on top of one other and separated by a 
non-absorbing nanoscale gap layer, such that the absorption spectra of each lattice when 
combined adds complimentarily (i.e., regions of low absorption in one lattice are compen- 
sated for by the high absorption of the other), the resonant-absorption characteristics of the 
PC augmented by partial disorder can potentially be exploited to the fullest extent possible 



for broadband absorption while also remaining feasible to large-scale industrial manufactur- 
ing. The low-index gap separation layer itself also provides additional mechanisms for light 
localization in nanostructured media that further contributes to enhancing absorption in the 
adjacent high-index layers while its effect on scattering-based Lambertian-textured films is 
marginal. A schematic of this design approach, somewhat exaggerated for illustrative pur- 
poses is shown in Fig. W[d) and (e) in which the tandem structure is first optimized for peak 
density in the spectrum (which amounts to maximizing the number of non-overlapping res- 
onant absorption modes of the two constituent PC lattices) and subsequently these narrow 
closely-spaced resonances are slightly broadened with the addition of disorder to create a 
more-uniform absorption profile that is large in amplitude, broadband and robust to incident 
radiation conditions. We consider here for simplicity the case of two PC slabs [4j in vacuum 
separated by an air gap which both incorporates all essential physical phenomena and has 
direct applications to conventional thin-film PV cells in which individual layers including 
both high-index semiconductors and low-index transparent conductive oxides are grown by 
thin-film deposition tools [251 [2S] . There is no need in the scope of the present work where 
the focus is solely on photon absorption in the silicon layer to include an anti-reflection (AR) 
coating in the front or a perfect reflector in the back as would be customary in an actual PV 
cell since the role of both components is mainly to enhance, oftentimes significantly [9l |27] , 
the absorption of existing resonances in the nanostructured films but not to give rise to new 
ones. Due to the complimentary way that the individual peaks of the two PCs combine in 
the tandem structure, the simplest lattice design consisting of a square lattice of circular 
holes etched into the silicon slab is adequate for good performance rendering unnecessary 
the need for intricate superlattice structures [131 [17] . Note however that in a real design the 
analysis of the overall absorption does not consist simply of a linear superposition of the two 
individual absorption spectra but rather is complicated by the various and more nuanced 
ways through which light propagates in the tandem structure (e.g., though some incident 
wavelengths pass through the top slab, they may be diffracted upwards by the bottom slab 
and potentially absorbed during a second pass) and also by the evanescent interaction be- 
tween the guided modes of the two slabs which is typically weak but non- negligible even 
when slightly separated as is the case here. 

To perform the topology optimization, we combine the capabilities of Meep, a freely- 
available open-source finite-difference time-domain (FDTD) tool [28], to compute the ab- 



sorption spectra at normal incidence with the nonhnear-optimization routines of NLopt [29] 
(details in Supplementary Information). Here we need not consider absorption at off-normal 
incidence since the addition of disorder in the final step will automatically reduce the sensi- 
tivity to incident radiation conditions [H] . Accurate topology optimization is made possible 
using Meep's subpixel averaging [201 ESI [30] which also significantly reduces the size of 
the computation by lowering the minimum spatial resolution required for reliable results. 
This is important since the objective function is iterated over a large number of times to 
explore small changes to geometrical parameters in order to engineer as many absorption 
resonances over the wideband spectrum as possible which tend to be highly susceptible to 
numerical artifacts introduced by "staircased" features jSHl EEl ED] enabling FDTD to be 
used as a versatile 3D design optimization tool. We use intrinsic crystalline silicon as our 
absorbing material and incorporate its full broadband complex refractive index profile [31] 
into the FDTD simulations with accuracy even near its indirect bandgap of llOSnm where 
absorption is almost negligible to obtain experimentally-realistic results (see Supplementary 
Information for more details). Although, for generality, modeling a tandem structure con- 
sisting of two completely-independent PC slabs with arbitrary lattice constants is preferred, 
incorporating two distinct unit cells into a single 3D simulation is computationally imprac- 
tical yet this is a minor design limitation as the other six structural parameters - as shown 
in Fig. [lie): the thickness of the top Si layer Vt (the bottom thickness V[, is known since the 
total thickness is fixed at l/im), the gap thickness g, the radius and height of the holes in the 
top and bottom lattices r^, ht, r^, hi, - provide sufficient flexibility for creating out-coupled 
Bloch-mode resonances over the entire range of the broadband solar spectrum (we also inves- 
tigated the computationally tractable case of two PC slabs with lattice constants differing 
by a factor of two though the results were not found to be an improvement to the single 
lattice-constant design). A planewave source is incident from above onto the tandem struc- 
ture and the absorption spectrum A{X), equivalent to 1-refiection-transmission, is calculated 
by Fourier-transforming the response to a short pulse. The absorptivity threshold used by 
the objective function to count the total number of peaks in the spectrum is taken to be 
that of our baseline performance metric: an equivalent lyum-thick Si film with Lambertian- 
textured surfaces [3l [32] which has an efficiency of 43.0% in our wavelength interval. Since 
resonances with especially-large peak amplitudes contribute most to increasing the overall 
efficiency when broadened, we add an extra 30% to our absorptivity peak threshold at each 



wavelength which while making the problem more challenging gives rise to better results. 
We impose no restrictions on the peak width or spacing relative to other peaks although 
these could potentially be used for further refinement. Once an optimal set of parameters for 
the two-lattice structure is determined by running multiple times with different randomly- 
chosen initial values to explore various local optima, we then form a supercell consisting of 
10x10 unit cells of the optimal structure and add positional disorder to each unperturbed 
hole in both lattices (while ensuring no overlap between holes to conserve the filling fraction) 
by an amount Api {Ap2) chosen randomly from a uniform distribution of values between 



and Api (Ap2) for both orthogonal in-plane directions for the top (bottom) slab. Three 
separate simulations are made for each structure and the results averaged due to the random 
nature of the design. 

Figure Wisu) shows the absorptivity spectra for three thin-film designs each with a total 
crystalline-silicon thickness of l/xm: an unpatterned film, a Lambertian-textured film [ob- 
tained from eq. (1) of Ref. [32] and finally the topology-optimized tandem design consisting 
of two PC slabs (top: thickness 708nm, hole radius 236nm, hole height 260nm and bottom: 
thickness 292nm, hole radius 199nm, hole height 244nm) with a lattice parameter of 641nm 
separated by a 228nm gap. The tandem design has numerous narrow, high-amplitude peaks 
- signatures of the coherent-resonant Bloch modes - that span the entire broadband spec- 
trum whereas the unpatterned slab has broad Fabry-Perot resonances with low amplitude. 
The complimentary way that the resonances of the individual slabs combine to form the 
tandem structure can be seen in Figure [2Fb) and (c) where the top slab accounts for most 
of the total number of peaks while the bottom slab contributes a few key resonances partic- 
ularly at longer wavelengths. Note that while the Lambertian-textured slab has little and 
diminishing absorption at long wavelengths where the absorption coefficient of crystalline 
silicon is vanishingly small the PC design, due to its resonant nature, has large absorption 
albeit appearing only as very-narrow peaks (due to the rate-matching phenomena discussed 
previously that underlies the resonant coupling between radiation and absorption by the 
material, as silicon's absorption coefficient becomes smaller at larger wavelengths leading 
to a corresponding increase in Tabs [3], the total lifetime of the resonant mode Ttot being 
1/Ttot=l/Tabs+1/Trad also iucrcascs resulting in the inversely-proportional peak width gener- 
ally becoming narrower which can be seen in Fig. [2]). Nevertheless, the tandem design with 
its maximized peak density, outperforms the Lambertian texture in light-trapping efficiency 



(48.8% versus 43.0%) although at off-normal angles and different polarizations of incident 
light the unperturbed lattices' resonance-based performance degrades significantly to below 
the Lambertian limit (both designs however are far better than the unpatterned film's 10.2% 
efficiency). Within the stacked arrangement, the efficiencies of the top and bottom slabs are 
42.6% and 6.2%, respectively, while in isolation they are 37.9% and 22.3% highlighting in 
part the importance of inter-slab interactions to the overall absorption of the tandem design. 
For comparison, the optimized lyum-thick single-slab design (lattice parameter, hole radius 
and height of 640nm, 256nm and 400nm, respectively) produces seven fewer resonances than 
the tandem design over the same wideband spectrum and therefore had a lower efficiency: 
nearly 5% below in absolute terms, yet still above the Lambertian limit. 

By proceeding to controUably introduce a partial amount of disorder into the topology- 
optimized tandem design, we can simultaneously boost efficiency and improve robustness to 
exceed the Lambertian limit by an even wider margin over a large angular cone. Fig. [3Fa) 
is a plot of the efficiency from eq. (pi) versus disorder for the optimized tandem-slabs and 
single-slab designs at normal incidence and shows that a positional disorder of approximately 



Api=Ap2=0.1a for the tandem slabs and Ap=0.15a for the single slab results in maximal 
light trapping of 9.8% and 6.6% above the Lambertian limit, respectively, while additional 
disorder beyond these partial amounts leads to a steady decrease of the efficiency in line with 
the analysis presented earlier. The tandem design therefore is roughly twice as effective as 
the single slab in overcoming the Lambertian limit due mostly to facilitating a larger number 
of absorption resonances. Note that because the two PC slabs that make up the tandem 
structure are different and the top slab in this case is the primary contributor to the overall 
absorption, the quantity of disorder applied to each lattice for enhancing light trapping need 
not be identical. We quantify the performance robustness of each design as the standard 
deviation of the efficiency averaged over normal (0°) incidence and five off- normal (10°, 
20°, 30°, 40°, 50°) angles of incidence for the TE polarization (electric field in the plane of 
incidence). Since more disorder results in better robustness [IH] which is a key requirement of 
a practical solar cell we make a slight trade off and apply not the quantities which maximize 



efficiency at normal incidence in the tandem design but values slightly greater {Api=0.2a 



and Ap2=0.25a) where the robustness is almost a factor of two larger: 50.7%±1.2% for the 
former versus 50.9%±0.7% for the latter. Fig. |3](b) demonstrates that this tandem design 
has greater absorption than the Lambertian scatterer at nearly every wavelength and angle 
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resulting in a light-trapping improvement of almost 10% above the Lambertian limit; the 
performance is also similar for the TM polarization. 

In summary, we have described a general design strategy derived from a new conceptual 
framework of photon capture for a nanostructured silicon thin film based on the quasi- 
resonant absorption of photons in a tandem arrangement of partially-disordered photonic- 
crystal slabs separated by a nanoscale gap where the overall light trapping surpasses a 
Lambertian-textured film by a wide margin over a large fraction of the solar spectrum and a 
broad angular cone. This demonstration of a photonic analogue to the mult i-j unction cell is 
made possible by recent advances in computational electrodynamics and may, with the full 
incorporation of electrical phenomena in the modeling as well as an experimental demon- 
stration, potentially lead to the development of high-efficiency and low-cost applications of 
renewable-energy systems for the harvesting of sunlight. 
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FIG. 1. (a) Complete random or so-called Lambertian textm'ing of the surface to isotropically 
scatter incident light rays into the plane of a weakly-absorbing film so as to increase the optical 
path length. First proposed nearly thirty years ago, this has thus far remained the most effective 
approach for light trapping over the widest set of conditions, (b) Photonic-crystal slab and other 
nanostructured designs in which light trapping occurs by resonant absorption into a guided mode 
depend on delicate wave-interference effects and are thus intrinsically narrowband and restricted to 
a small angular cone. Strong disorder improves robustness but nevertheless comes at the expense of 
overall light trapping relative to the unperturbed case and is still less effective than the Lambertian 
texture, (c) Tandem arrangement of two PC slabs, both consisting of a square lattice of circular 
air holes in silicon, stacked vertically on top of each other. Shown are the degrees of freedom - slab 
thicknesses vi and V2, lattice parameters ai and 02, hole radii ri and r2, hole heights hi and h2 
and gap separation g - used in the topology optimization to (d) engineer as many non-overlapping 
resonances over the wideband solar spectrum as possible. Following this, (e) each hole is perturbed 
from its position in the unperturbed lattice by amounts Api (Ap2) chosen randomly from a uniform 



distribution of values between and Api (Ap2) for both orthogonal in-plane directions of the top 
(bottom) slab to boost light trapping and robustness by creating a more- uniform absorption profile. 
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FIG. 2. (a) Absorption versus wavelength profile at normal incidence for three thin- film photo- 
voltaic designs each with a total intrinsic-crystalline-silicon thickness of 1/im: an unpatterned slab 
(green), a Lambertian-textured slab (red) [obtained from eq. (1) of Ref . [32] and the topology- 
optimized tandem PC slabs (blue). The tandem PC slabs both consist of a square lattice (pe- 
riodicity, a=641nm) of holes in silicon separated by a 228nm gap. Shown for each design is the 
photon-absorption efficiency defined in eq. ([3| which is a measure of light trapping relative to 
a perfect absorber, (b) and (c) Individual absorption spectra for the top (slab #1: wt=708nm, 
rt=236nm, /it=260nm) and bottom (slab #2: Wb=292nm, rfe=199nm, /i{,=244nm) PC slabs of the 
optimized tandem design demonstrating how the resonances add compliment arily over the broad- 
band spectrum. 
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FIG. 3. (a) Light-trapping efficiency from eq. ([s]) as computed from the absorption profile at 
normal incidence versus hole positional disorder for the topology-optimized tandem slabs (blue) 



and single slab (red) showing that partial disorder (tandem: Api=Ap2=0-la, single: Ap=0.15a) 
maximizes the light trapping (tandem: 52.8%, single: 49.6%) while additional disorder is sub 
optimal and leads to a steady decline. Note that the tandem design is nearly twice as effective as 
the single slab design in surpassing the Lambertian limit, (b) Absorption versus wavelength profile 
at normal (0°) incidence and five off-normal (10°, 20°, 30°, 40°, 50°) angles of incidence for the 
TE polarization (electric field in the plane of incidence) of the optimized tandem design with the 



addition of partial disorder of Api=0.2a and Ap2=0.25a as shown in the inset schematics. The 
absorption profile for the individual angles are colored while the averaged profile is shown in black. 
The absorption of the optimized tandem PC slabs with partial disorder exceeds the Lambertian 
texture for nearly every wavelength and angle resulting in an overall light-trapping efficiency that 
is approximately 10% greater. The TM polarization data (not shown) has similar characteristics. 
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